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Abstract. This talk 1 presents work concepts and results for the determination of the 
fine structure constant a at the Zq mass resonance. The problem consisting of the 
break-down of global duality for singular integral weights is circumvented by using 
a polynomial fit which mimics this weight function. This method is conservative in 
the sense that it is mostly independent of special assumptions. In this context the 
difference between local and global duality is explained. 

INTRODUCTION 

There is a great deal of interest in the accurate determination of the running fine 
structure constant a at the scale of the Z$ mass [1-4]. The value of oc(Mz) is of 
paramount importance for all precision tests of the Standard Model. Furthermore, 
an accurate knowledge of aiMz) is instrumental in narrowing down the the mass 
window for the last missing particle of the Standard Model, the Higgs particle. 

The main source of uncertainty in the determination of a(Mz) is the hadronic 
contribution to the e + e~ annihilations needed for this evaluation. The necessary 
dispersion integral that enters this calculation has in the past been evaluated by 
using experimental e + e~ annihilation data. Discrepancies in the experimental data 
between different experiments suggest large systematic uncertainties in each of the 
experiments. In order to reduce the influence of the systematic uncertainties on 
the determination of a(Mz) one may attempt to add some theoretical input to the 
evaluation of the hadronic contribution to a(Mz)- 

M. Davier and A. Hocker [3] use QCD perturbation theory in form of local duality 
(explained later on) in the region above s = (1.8 GeV) 2 for the light flavours, while 
J.H. Kiihn and M. Steinhauser [5] use perturbative results for energy regions outside 
the charm and bottom threshold regions. Our approach [6] is quite different. We 
attempt to minimize the influence of data in the dispersion integral over the whole 
energy region including the threshold regions. 
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Local and global duality 



For the e + e~ annihilation process there is a main connection between the spectral 
density p(s) and the two-point correlator IT(g 2 ) given in kind of the dispersion 
relation 
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which implies the reverse relation 

p(s) = ^-DiscTl(s) (2) 

where the discontinuity is given by 

Disc U(q 2 ) = Tl(q 2 e- i7T ) - U(q 2 e i7T ) (3) 

and So — 4m 2 is the production threshold of the light flavours. These relations 
can be a chain between the theory, i.e. the two-point correlator function within 
perturbative QCD on the one hand and the experiment, i.e. the spectral density or, 
equivalently, the total cross section on the other hand. But there is one obstacle in 
using these relations: As depending on methods of functional analysis, the inverse of 
the dispersion relation is only valid if there are no poles encircled by the path in the 
complex plane, a condition which is necessary to obtain this relation. These poles 
can have their origin from weight functions in combination with the spectral density. 
This means that if there is such a weight function included in the integration of 
the spectral density, the inverse relation shown above is only valid locally and not 
globally. We call this local resp. global duality. 

The experiment side 

The hadronic contribution to a(Mz) which we are concentrating on is given by 
the integral [1] 

a f°° 

A« had (M z ) = —Re / R(s)H(s)ds (4) 
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where R(s) is the total e + e~ hadronic cross section and H(s) is the weight function 

M 2 

The hadronic cross section is related to the spectral density by 

R(s) = 12tt 2 p(s). (6) 

But we see: the weight function H(s) is indeed singular at the points s = and 
s = Mf on the real axis. So global duality is not valid any more. 



The theory side 



The two-point correlator is given by 

i J (OUTW^Me^x = (-g af sq 2 + q a q P )U(q 2 ) 



(7) 



where we only included the isospin contribution I — 1, in contrast to corresponding 
considerations for the r decay. The scalar correlator function II(g 2 ) consists of 
perturbative and non-perturbative contributions which we include to the extend we 
need them to keep the accuracy. For the perturbative contribution to the correlator 
we use a result given in ref. [7]. I only write down the first few terms, 
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with L = ln(/x 2 /g 2 ) while in ref. [7] the expression is given up to 0(a 2 ,m 1 q 2 /q 12 ). 
The number of active flavours is denoted by rif. For the zeroth order term in the 
m 2 /q 2 expansion we have added higher order terms in a s , 
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with k — 1, ki — 1.63982 and k 2 = 6.37101. We have denoted the yet unknown 
constant term in the four-loop contribution by C3. Remark, however, that the 
constant non-logarithmic terms will not contribute to our calculations. The non- 
perturbative contributions are given in ref. [8], 
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where we have included the m 4 /g 4 -contribution arising from the unit operator. In 
this expression we used the SU(3) colour factors Cf = 4/3, Ca = 3, and T F = 1/2. 
For the coupling constant a s as well as for the running quark mass we use four-loop 
expression given in refs. [9-11] even though in both cases the three-loop accuracy 
would already have been sufficient for the present application. 



INTRODUCING OUR METHOD 



Our method is based on the fact that we can use global duality when the weight 
function is non-singular. This is the case for a polynomial function. So we mimic 
the weight function by a polynomial function obeying different conditions which 
we will explain later. By adding and subtracting this polynomial function Pn(s) 
of given order N to the weight function H(s), we obtain without any restrictions 

f Sb p(s)H(s)ds = H p(s) (H(s) - P N (s)) ds + H p(s)P N (s)ds (11) 

J s a J S a J s a 

where [s a , Sf,] is any interval out of the total integration range. But because the 
second term has now a polynomial weight, we can use global duality to write 

P^b ^ P^b 

/ p{s)P N {s)ds = — ;/ DiscIl(s)P N {s)ds = 

Js a 2%i Js a 

= -7^$ l n(-s)P N (s)ds + ^-£ U(-s)P N (s)ds. (12) 

2m J\s\=s a 2m J\s\=s b 

Thus this part can be represented by a difference of two circle integrals in the 
complex plane. On the other hand, the difference H(s) — Pn(s) suppresses the 
contribution of the first part. Our method consists thus of the following steps: 

• replacing p(s) in the first part of Eq. (11) by the value of the experimentally 
measured total cross section R(s) (see e.g. ref. [1]) 

• replacing the circle integral contribution to flavours at their threshold by zero 

• in all other cases inserting the QCD perturbative and non-perturbative parts 
of II(— s) on the circle 

These replacements can be seen as a concept within QCD sum rules. To obtain the 
best efficience of our method, we have to restrict the polynomial function by the 
following contraints: 

• The method of least squares should be used to mimic the weight 

• However, the degree N should not be higher than the order of the highest 
perturbative resp. non-perturbative contribution increased by one (this is a 
consequence of the Cauchy's theorem which is involved in the analytical inte- 
gration of the circle integrals) 

• Especially for the low energy region, the polynomial function should vanish 
on the real axis to avoid instanton effects 

• In regions where resonances occur, the polynomial function should fit the 
weight function to suppress those contributions which constitute the highest 
uncertainty of the experimental data 



As just mentioned, the integration on the circle can be done analytically by using 
the Cauchy's theorem. But we have to keep in mind that the result for II(— s) 
we use here depends logarithmically on the renormalization scale \x and on the 
parameters of the theory that are renormalized at the scale \x. These are the 
strong coupling constant, the quark masses and the condensates. As advocated 
in [12], we implement the renormalization group improvement for the moments of 
the electromagnetic correlator by performing the integrations over the circle with 
radius \s\ = Sb with constant parameters, i.e. they are renormalized at a fixed scale 
/i. Subsequently these parameters are evolved from this scale to fi 2 = s a using the 
four-loop P function. In other words, we impose the renormalization group equation 
on the moments rather than on the correlator itself. This procedure is not only 
technically simpler but also avoids possible inconsistencies inherent to the usual 
approach where one applies the renormalization group to the correlator, expands 
in powers of ln(s//x 2 ) and carries out the integration in the complex plane only at 
the end. In the present case the reference scale is given by A^g. 



Subdividing the integration interval 

As a first interval we select the range from the light flavour production threshold 
s = 4m 2 and the next threshold marked by the mass of the ip, s± — m 2 ~ 
(3.1 GeV) 2 . In this case we set the inner circle integral to zero and obtain 

f 1 R(s)H(s)ds = f * R exp (s) {H{s) - P N (s)) ds + 6vri <f U QCD {-s)P N {s)ds. 

J so J so J\s\=si 

(13) 

As mentioned above, we impose the constraints to the polynomial function that it 
should vanish on the real axis at s — si and should coincide with the weight function 
at the p resonance, i.e. for s = m 2 ~ (1 GeV) 2 . Fig. 1 shows polynomials of different 
order in comparison with the weight function. The results shown in Fig. 2 are 
compared with the result obtained by using only the experimental data. For the up 
and down quarks we only keep the mass independent part of the QCD contribution 
while for the strange quark we include also the terms to order 0(m 2 s /q 2 ). 

The second interval is limited by Si and the threshold marked by the mass of 
the T, S2 = m\ ~ (9.46 GeV) 2 . For the charm quark, we again set the inner circle 
integral to zero, but for the lighter quarks we have to keep both. The perturbative 
series for the charm quark is used up to it's known extend. 

The third interval given between S2 and (40 GeV) 2 is again subdivided into 
two pieces because of it's length. For the first of these two intervals we choose 
[(9.46 GeV) 2 , (30 GeV) 2 ], for the second [(30 GeV) 2 , (40 GeV) 2 ]. Now the bottom 
quark is the one for which the "threshold rule" (i.e. leaving out the inner circle) 
applies. The remaining part of the integral starting from s 4 = (40 GeV) 2 up to 
infinity is done using local duality, i.e. by inserting the function R(s) obtained for 
perturbative QCD into the second part of Eq. (11). 
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FIGURE 1. Weight function H(s) and polynomial approximations Pn(s) in the lowest energy 
interval Im^ < -Js < 3.1 GeV. The least square fit was done in the interval m p < y/s < 3.1 GeV 
with further constraints H(s) = Pn{s) at ^/s = 1 GeV and Pn(s) = at ^/s = 3.1 GeV. The qual- 
ity of the polynomial approximations are shown up to N — 4. We use the scaled variable s/si for 
the polynomial approximation where si is the upper radius such that P/v(s/si) is dimcnsionless. 
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FIGURE 2. Comparison of the l.h.s. and r.h.s. of the sum rule given by Eq. (13) in the interval 
0.28 GeV < y/s < 3.1 GeV. Dotted horizontal line: value of integrating the l.h.s. using experimen- 
tal data including error bars [1]. The points give the values of the r.h.s. integration for various 
orders TV of the polynomial approximation. Straight line interpolations between the points are 
for illustration only. The dashed lines indicate the error estimate of our calculation. 



Our results are collected in Table I. To obtain these results, we used the con- 
densate values 

(-GG) = (0.04 ± 0.04) GeV 4 , a s (qq) 2 = (4 ± 4) • 10~ 4 GeV 6 . (14) 

7T 

For the errors coming from the uncertainty of the QCD scale we take 

A m = 380 ± 60 MeV (15) 

The errors resulting from the uncertainty in the QCD scale in different energy 
intervals are clearly correlated and will have to be added linearly in the end. We 
also include the error of the strange quark mass in the light quark region which is 
taken as 

m s (lGeV) = 200 ±60 MeV (16) 

For the charm and bottom quark masses we use the values 

fh c (m c ) = 1.4 ± 0.2 GeV, rh b (m b ) = 4.8 ± 0.3 GeV. (17) 

Summing up the contributions from the five flavours u, d, s, c, and b, our result for 
the hadronic contribution to the dispersion integral including the systematic error 
due to the dependence on (column 5 in Table 1) reads 

Aa^ d (M z ) = (277.6 ± 4.1) • 10~ 4 . (18) 

In order to obtain the total result for a(M z ), we have to add the lepton and 
top contributions. Since we have nothing new to add to the calculation of these 
contributions we simply take the values from ref. [5], 

Aa Ld( M z) = (-0.70 ± 0.05) • 10^ 4 , Aa lcp (M z ) « 314.97 • 10~ 4 . (19) 

Writing Aa(M z ) = Aa lcp (M z ) + Aa had (M z ) our final result is (a(O)" 1 = 137.036) 

a(M z )~ 1 = a(0) _1 (l - Aa(M z )) = 128.925 ± 0.056. (20) 



CONCLUSION AND OUTLOOK 

I have presented a method to obtain the running fine structure constant a at the 
scale of the Z mass with minimal input of experimental data. This method is con- 
servative in the meaning that its error is as free from assumptions as possible. Our 
method is discussed and compared with other methods (see e.g. ref. [13] - however, 
with our preliminary results). In ref. [14] Matthias Steinhauser says that "it is very 
impressive that the new analysis show very good agreement both in their central 
values and their quoted errors." I nevertheless would like to close this talk with 
the remark that all recent calculations of a(M z ) should not deter experimentalists 
from remeasuring the e + e~ annihilation cross section more accurately in the low 
and intermediate energy region, as such data are absolutely essential for a precise 
value of a(M z ), unbiased by theory. 



TABLE 1. Contributions of different energy intervals to a^ d (M z ). Second col- 
umn: choice of neighbouring pairs of the polynomial degree N. Third column: 
fraction of the contribution of experimental data [1]. Fourth column: contribution 
to Aa^ d (Mz) with all errors included except for the systematic error due to the 



dependence on A^g- 


which is 


separately listed 


in the fifth column. 






interval 


values 


data 


contribution 




error 


for i/s 


of N 


contribution 


to Aa^ d (M z ) 


due to Aj^g 


[0.28 GeV, 3.1 GeV] 


1,2 


24% 


(73.9 ± 1.1) • 10" 4 


0.9 


10~ 4 


[3.1 GcV, 9.46 GeV] 


3,4 


0.3% 


(69.5 ±3.0) • 10~ 4 


1.4 


io~ 4 


[9.46 GeV, 30 GeV] 


3,4 


1.1% 


(71.6 ±0.5) • 10~ 4 


0.06 


10~ 4 


[30 GcV, 40 GcV] 


3,4 


0.15% 


(19.93 ±0.01) • 10~ 4 


0.02 


10~ 4 


y/s > 40 GcV 






(42.67 ±0.09) • 10~ 4 






total range 






(277.6 ±3.2) • 10~ 4 


1.67 


10~ 4 
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